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Abstract 

Models of inviscid incompressible fluid are considered, with the kinetic energy (i.e., 
the Lagrangian functional) taking the form ~ / /c^lv^pd^k in 3D Fourier represen- 
tation, where a is a constant, < a < 1. Unlike the case a = (the usual Eulerian 
hydrodynamics), a flnite value of a results in a finite energy for a singular, frozen- 
in vortex filament. This property allows us to study the dynamics of such filaments 
without the necessity of a regularization procedure for short length scales. The linear 
analysis of small symmetrical deviations from a stationary solution is performed for a 
pair of anti-parallel vortex filaments and an analog of the Crow instability is found at 
small wave-numbers. A local approximate Hamiltonian is obtained for the nonlinear 
long-scale dynamics of this system. Self-similar solutions of the corresponding equa- 
tions are found analytically. They describe the formation of a finite time singularity, 
with all length scales decreasing like {t* — tY^^'^~°'\ where t* is the singularity time. 

PACS numbers: 47.15.Ki, 47.32.Cc 

1 Introduction 

The question of the possibility for the spontaneous formation of a finite time singularity 
in solutions of the Euler equation for an ideal incompressible fluid has been discussed for 
a long time. At present, this fundamental problem of fluid dynamics is still far from a 
complete solution, though some rigorous analytical results have been obtained. So, it has 
been found by Beale, Kato and Majda that no singularity occurs if the time- integral of the 
maximum of the vorticity magnitude is finite. Another result (Constantin and co-workers, 
0"[S]) postulates that a blow-up of the vorticity, if does take place, must be accompanied 
by a singularity in the field of the vorticity direction. In general, the nature of the presumed 
singularity has not yet been clarified, although many theoretical scenarios for blow-up have 
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been suggested until now, and also extensive numerical simulations have been performed to 



observe the singular behavior (see 0]-[|T^ and references therein). In particular, the locally 
self-similar regime of singularity formation seems very probable. For this regime, a region of 
finite size may be distinguished in flow, where all length scales, corresponding to vorticity 
distribution, decrease like {t* —tY^"^, the velocity increases according to {t* — and the 
maximum of the vorticity behaves like {t* — t)^^. It is necessary to emphasize that this is the 
only possible scaling, which is compatible both with the dimensional structure of the Euler 
equation and with the freezing-in property of the vorticity. In this process, an accelerated 
straining of vortex lines takes place, and it is the reason for amplification of the vorticity 
magnitude. It is very important that the curvature of vortex lines in the assumed self-similar 
solutions should tend to infinity in the vicinity of the singular point, in accordance with the 
result of Constantin, Fefferman and Majda. Thus, the problem of singularity in the curvature 
of vortex lines and the problem of singularity in their local stretching are closely connected. 

In this paper, we take the point of view that infinite curvature of frozen-in vortex lines 
is in some sense a more fundamental characteristics of hydrodynamic singularity than in- 
finite value of the vorticity maximum. To illustrate this statement, we consider a class of 
models of an incompressible inviscid fluid, different from Eulerian hydrodynamics, such that 
finite energy solutions with infinitely thin frozen-in vortex filaments of finite strengths are 
possible. Thus, we deal with a situation when the vorticity maximum is infinite from the 
very beginning, but nevertheless, this fact itself does not imply a singular behaviour in the 
dynamics of vortex strings, while their shape is smooth and the distance between them is 
finite. However, the interaction between filaments may result in formation of a finite time 
singularity for the curvature of vortex strings. It is the main purpose of present work to 
study this phenomenon analytically. 

It is a well known fact that absence of solutions with singular vortex filaments in Eulerian 
hydrodynamics is manifested, in particular, as a logarithmic divergency of the corresponding 
formal expression for the energy functional of an infinitely thin vortex filament having a finite 
circulation F and a shape R(0 (this is actually the Hamiltonian functional determining 
entirely the dynamics of the system, as is shown in Section 2): 

«W«» = 5^// |R(6)-R(&)| 

More important is that the self-induced velocity of a curved string in Eulerian hydrodynamics 
is also infinite. This is the reason, why we cannot work in the framework of Eulerian 
hydrodynamics with such one-dimensional objects, that are very attractive for theoretical 
treatment. The situation becomes more favourable, when we consider a class of regularized 
models, with the divergency of the energy functional eliminated. It should be stressed here 
that in regularized systems the usual relation ft = curl v between the vorticity and velocity 
fields is no more valid, and in this case F is not the circulation of the velocity around 
the filament, but it is the circulation of the canonical momentum field (see Section 2 for 
more details). However, dynamical properties of a de-singularized system depend on the 
manner of regularization. For instance, it is possible to replace the singular Green's function 
G(|Ri — R2I) in (lU) (where G{r) ~ 1/r) by some analytical function which has no singular 
points near the real axis in the complex plane (for examples by Gq{r) ~ tanh(gr)/r or by 
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G^{r) ~ l/Vr^ + e^). In that case we may not expect any finite time singularity formation, 
because tlie corresponding velocity field created by the vortex string appears to be too 
smooth with any shape of the curve, and this fact prevents drawing together some pieces of 
the string. With such a very smooth velocity field, a singularity formation needs an infinite 
time. 

In this paper we consider another type of regularization of the Hamiltonian functional, 
when the Green's function is still singular, but this singularity is integrable in the contour 
integral analogous to the expression (|lD: 

-wrjRmi // (RXei)-R^(6))cieic^6 ... 
HJR(0} -Yff |R(^0-R(6)|-° ' ^ ^ 

with a small but finite positive constant < a <^ 1. If a is not small, we actually have 
models that are rather different from Eulerian hydrodynamics. Nevertheless, such mod- 
els still have many common features with usual hydrodynamics, which are important for 
singularity formation in the process of the interaction between vortex filaments: a similar 
hydrodynamic type structure of the Hamiltonian and a power-like behaviour of the Green's 
function, with negative exponent. Therefore we believe that it is useful to investigate these 
models, especially the question about the formation of a finite time singularity in the vortex 
line curvature. We hope the results of our study will shed more light on the problem of 
blow-up in Eulerian hydrodynamics. 

This paper is organized as follows. In Section II, we briefly review some basic proper- 
ties of frozen-in vorticity dynamics in a perfect fluid, with giving necessary definitions for 
theoretical conceptions used in our study. In general, our approach is based on the Hamil- 



tonian formalism for frozen-in vortex lines |T5-|T^. Then, in Section HI, we perform the 



linear analysis of stability for a pair of symmetric anti-parallel vortex filaments and find an 
instability at small wave numbers analogous to the Crow instability [T^. In Section IV, 
we postulate a local approximate Hamiltonian for the long scale nonlinear dynamics of the 
pair of filaments and present analytical self-similar solutions of the corresponding equations. 
Those solutions describe finite time singularity formation, with the length scales decreasing 
like (t* - |-}^ig ig i-j^g 

main result of present work. In Section V, we make some 
concluding remarks about vortex filaments of a finite width, about long scale approximation 
for systems with the Green's function of a general form, and finally about how it is possible 
to improve the approximation in the case of small a, when the unstable region is narrow in 
wave number space. In Appendix A, we write in terms of the special mathematical functions 
some integral expressions needed for calculation of instability increment of the vortex pair. 
In Appendix B, we provide details about the integration procedure for the system of ordinary 
differential equations related to the self-similar solutions. 



2 Hamiltonian dynamics of vortex filaments 

To clarify the meaning of the suggested models @ and to explain the employed theoretical 
method, we recall some general properties of frozen-in vorticity dynamics in a perfect fluid. 



starting from the Lagrangian formalism 1T^-||2^, [|T4|-[jl7 . 
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Let a Lagrangian functional £{v} specify the dynamics of some incompressible medium 
of unit density, with the solenoidal velocity field v(r, t). Especially we are interested here in 
systems with quadratic Lagrangians, which in 3D Fourier representation take the form: 



'Cm{'v} 



_fk_ 

(2^ 



M(A;)|vk| 



(3) 



where M{k) is some given positive function of the absolute value of the wave vector k. This 
expression should be understood as the kinetic energy on the group of volume preserving 
mappings x(a, t), and the velocity field v(x, t) is defined as the time derivative x(a, t) taken 
at the point a(x, t). Obviously, all the systems (Q) possess the properties of homogeneity and 
isothropy in the space. It is clear that the usual Eulerian hydrodynamics corresponds to the 
simplest case M{k) = 1. Another physically important example concerns the homogeneous 
incompressible electron magnetohydrodynamics (EMHD), for which M{k) = 1 + with 
a constant q beeng the screening parameter [|l^,||26|. Also the case M{k) = 1 + A^A;^ 
has been studied, with a constant A, which corresponds to the so called avereged Eulerian 
hydrodynamics (see, for instance, the papers ||2^, for more details). In the general 
case, the systems (§) may be understood as models for some inviscid non-newtonian fluids. 
It should be noted that there exists a direct relation between such models and the vortex 
blob method introduced by Chorin for desingularization of the Eulerian hydrodynamics [^ |. 
Some discussion of this relation, for the case of the averaged Eulerian hydrodynamics, can 
be found in the papers p^] and |^ 



Due to the presence of the Noether type symmetry with respect to relabeling of La- 
grangian labels of fluid points [^-|^, |1^-|T6[, all such systems have an infinite number 



of integrals of motion, which can be expressed as conservation of the circulations of the 
canonical momentum field p(r,t), 

along any closed contour c(t) advected by flow, thus the generalized theorem of Kelvin is 
valid: 

Tc = 9 (p ■ d\) = const. (5) 

Jc{t) 

These integrals of motion correspond to the frozen-in property of the canonical vorticity field 

O = curl p = curl — . (6) 
ov 

After defining the Hamiltonian functional 7Y{f2}, 
the equation of motion for the vorticity takes the form 



v=v{n}' 



(7) 



fit = curl 



curl 




X O 
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This equation describes the transport of frozen-in vortex hnes by the flow having the velocity 
field 

v = curl(^g), (9) 

It is very important in this process that all topological characteristics of the vorticity field 
are conserved |l^, It follows from Eqs. © and that the Hamiltonian TIm 

corresponding to the Lag rangian JUm is 

with the Green's function Gm{''^) being equal to the following integral: 

^ , , , rf^k e'^"" 1 ^r'sinfcr dk 



(27r)3 A;2M(A;) 27r2 7 kr M{k)' 

The frozen-in vorticity field can be represented in topologically simple contin- 
uous distribution of vortex lines [0]-[|l^: 

n{r, t)= rfV j 5{v - R(z/, t))^di, (12) 

where a 2D Lagrangian coordinate v = (z^i, ^2)1 which lies in some manifold A/", is the label 
of a vortex line, while the longitudinal coordinate ^ determines a point on the line. 

The important characteristics of the system: the (virtual) linear momentum P and the 
angular momentum M can be expressed as follows: 



P 

M 



J^d^iy^f[RxR^]d^, (13) 
J^d'u y[Rx[RxR^]]dC (14) 



In the limit, when the shapes R(z/, ^,t) of vortex lines do not depend on the label z/, we 
have one singular vortex filament with a finite circulation F = Jj^d'^u. In this case, the flow 
is potential in the space around the filament: p = V$, with a multi-valued scalar potential 
^(r, t). The potential flow domain is passive from the dynamical viewpoint, because there 
the flow depends entirely on the filament shape. The dynamics of the shape R(^, t) of such 
infinitely thin vortex filament is determined in a self-consistent manner by the variational 
principle with the Lagrangian £^j{R} p!4|-||r7|, 

/:^ = r^([R'xRi].D(R))rfe-y^^GA/(|R(6)-R(6)l) (R'(6) ■ R'(6)) c^«2, (15) 

where the vector function D(R) must have unit divergence |T^: 

divRD(R) = 1. (16) 
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The generalization of the expression (15) to a case of several filaments with the circulations 
r^") and shapes Il^^\^,t), n = 1..N, is straightforward: one should write a single sum over 
n for the first term and a double sum for the Hamiltonian. 

It is easy to see that the Hamiltonian (|^) corresponds to the function M{k) in the form 

M{k) ~ A;". (17) 

The choice of the longitudinal parameter ^ is not unique, but this does not affect the dy- 
namics of the vortex string which is an invariant geometric object. Sometimes it is convenient 
to use parameterization of the vortex line shape by the Cartesian coordinate: 

R(e,t) = (x(e,t),F(e,t),o- (is) 

Then the choice D = (0, F, 0) gives immediately that X(^,t) and Y{C,,t) are canonically 
conjugated quantities. 

Hereafter, we will consider vortex filaments with unit circulation for simplicity. So 
the symbol F, if appearing in some expressions below, will mean the special mathemati- 
cal Gamma function. Also, without loss of generality, all quantities may be considered as 
dimensionless. 

Now, for some fixed value of the parameter a, let us consider the symmetrical dynamics 
of a pair of oppositely rotating vortex filaments, with the symmetry plane y = const. Due to 
this symmetry, it is sufficient to consider only one of the filaments. It follows from the above 
discussion that the exact expression for the Hamiltonian of this system is the following: 

Iff {l + X[X^ + YIY{)d^,d^2 

rir 



2 



((^i-6)^+(^i-^2)2+(ri-F2)2) 

1 ff (-i-xjx^ + y/y^') di^di^ 



l-a 
2 



2 



((^1 - 6)2 +(^1 - X^Y +(Fi + F2 + hf 



l — a 1 
2 



(19) 



where h is the mean distance between the two filaments {h does not depend on time because 
of the conservation law for the momentum (|13])), Xi = X(^i), X[ = X'(^i) and so on. The 
first term in Eg . ([19|) describes the non-local self-interaction of the filament, while the second 
one corresponds to the interaction with the second filament. The Hamiltonian equations of 
motion have the form 

y(0--§^y (20) 



3 Crow instability for a pair of vortex filaments 

The system with the Hamiltonian (|1^) possesses the exact stationary solution 

X(e,t) = C(a,6)t, r(^,t) = 0, (21) 

which describes the uniform motion of straight filaments. Here the stationary velocity C(q;, h) 
is proportional to But this solution appears to be unstable due to an analog of the 
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Crow instability |18]. In this section we consider the hnear evolution of small perturbations 
of the stationary solution, and derive the linear growth rate. 

To perform the linear analysis of small deviations of the vortex shape from a straight 
line, we need the quadratic part of the Hamiltonian ([191) : 



-2JJ + 2 J J [^) 16-6'^-" 



1 /•/• /«-3\ 2h'^{Yi + Y2f 



2 



5-a 
2 



_diidi2. (22) 

;(ei - ^2? + h^] 

For further consideration, it is useful to rewrite it in the ID Fourier representation: 

n^I^ = \j ^(A^{k)X,X_, + B^{k)Y,Y_,). (23) 

Expressions for the functions Aa{k) and Ba{k) follow from Eq.(p^. So, Aa{k) can be 
represented as follows: 

f+oo / 1 1 \ 

+2(a-l)r-y^ (1 - cos(A:6C)) [-^-j^^^^j 

= 2(1-«)V-M (l-cos(fc6C)) 3^^+ ^ sr^ MC (24) 

Obviously, Aq,(/c) is positive everywhere. Analogous calculations for the function Ba{k) give: 

dC 



r+oo ( \ \ \ H 

BM = 2k%'^ I cosikbC) + J^^^Y)^) dC + 2(a - 1)&"-^ 



:i-cos(A;6C)) 



3-a 



+2(1 _ + eos(**C)) f^^j^ + ^^^) iC 

dC 







= 2(1 - afb'^-' ll - cos(fc6C))^ 

-2(l-«)(3-«)6"-2/ (l + cos(fc6C)) ( W-- - ^25) 

In Appendix A, Aq,(A;) and Ba{k) are expressed through the Euler Gamma function T{x) 
and the modified Bessel functions of the second kind K^ix). 
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Figure 1: The dependences u:l{k) = Aa{k)Ba{k) with 6 = 1 for a = 0.01, 0.025, 0.05, 0.1, 0.25, 
0.5. Lines corresponding to the given values of a intersect the horizontal axis in the indicated order. 



The dispersion relation between the frequency oOa of a small amplitude perturbation of 
the filament shape and the corresponding wave number k is simply given by the formula 

ulik) = A^{k)B^{k), (26) 

since the linearized equations of motion for and are 

Xk = B^{k)Yk, n = -A^{k)Xk, (27) 

as follows from Eq. (pO|) . In Fig.|I]we have plotted tu^ versus k for several values of a. 

It is easy to see that at small wave numbers the product Aa{k)Ba{k) is negative. Indeed, 
after some calculations we obtain in leading order for kb 1: 

AM ^ {^^) h-a, (28) 



a 



5„(fc) ^ -4(1 - aYb'^-'h-a, (29) 
where the constant I^-a is given by the integral 



'3-a 



iC + i)"^ 2r(^) ' 



with r(..) being the Gamma function. Therefore, an instability takes place at small k. The 
unstable domain in the wave number space corresponds to a range < qoi<y) where Ba{k) 
is negative, with the function qo{a) behaving, at small values of a, like -y/a: 

qo{a) ^ 2y/a, a < 1. (31) 



The plot of qo{a) is shown in Fig.|[ The instability increment 7a(fc) = J —Aa{k)Ba{k) is 
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Figure 2: The boundary of instability qo{a). 

proportional to the absolute value of k at very small values of kb: 

7(A;) ^ (1 - a)h-a ■ 2|A;|6"-V(1 - (32) 

However, for each a there exists a maximum value 7max(ct) of the increment, which is attained 
at kb ~ ^/a. Therefore the approximate expressions (|29|) and (|3^) are valid only if | 1 6 <^ ^/a. 

For large wave numbers, \k\b ^ 1, the functions Aa{k) and Ba{k) are both positive. The 
asymptotic approximations in that region are: 

Mk) « BJk) « 2(1 - r ^i^rf, = ^._„ 2(l-a)cos(W2)r(a) ^ 

JO r]'^ 2 — a 

Note that this expression does not contain the parameter b. For a single vortex filament it 
is actually the exact expression for Aa{k) and Ba{k), which is valid in the whole range of k. 

A general nonlinear analysis of the non-local system (|19D is difficult. Therefore we need 
some simplified model which would approximate the nonlinear dynamics, at least in the most 
interesting long scale unstable regime. In the next section, we suggest such an approximate 
model and find a class of solutions describing the formation of a finite time singularity. 



4 Singularity in long-scale nonlinear dynamics 



We note that the same long-scale limit as (p^-p9|) can be obtained from the local nonlinear 
Hamiltonian 

{2YY 



HiiRiO} = (1 - «)/3-« f + Z'^d^, (34) 

J a 

where the coordinate Y{^) is measured from the symmetry plane. This Hamiltonian approx- 
imates the exact non-local Hamiltonian of a symmetrical pair of vortex filaments in the case 
when the ratio of a typical value of F to a typical longitudinal scale L is much smaller than 
qo{a): 

Y/L < v^. (35) 
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In particular, this means that the slope of the curve with respect to the symmetry plane 
should be small, and also Y should be small in comparison with the radius of the line 
curvature. When Y = const, X' = const, Z' = const, expression ( PH ) gives the same result 
for uniform stationary motion as the exact Hamiltonian. 

With the Cartesian parametrization (|1^), the corresponding approximate local nonlinear 
equations of motion have the form (after appropriate time rescaling) 

• 1 VTTx^ 

^ = (2 -a) ' 

^ = J2^{^l + X'^) ^^^^ 

and they allow to obtain a simple explanation of the instability. On a qualitative level 
of understanding, the reason for the instability is that if initially some pieces of the curve 
were more close to the symmetry plane and convex in the direction of motion, then at 
subsequent moments in time the curvature will be increasing because of smaller values of Y 
and corresponding larger velocity, while Y will be decreasing due to the curvature. Thus, 
the feedback is positive and the system is unstable. In the final stage of the instability 
development, a locally self-similar regime in the dynamics is possible, because the above 
equations admit the self-similar substitution 

X{^,t) = X*-{f-tfx{i^-C){t*~t)-^), (38) 
F(e,t) = it*-try{i^~C)it*-tr^), (39) 

with arbitrary constants X*, if:*, and with the exponent 



After substituting Eqs. ( |55| - PPD into Eqs.( PB[ - P7| ), we obtain a pair of ordinary differential 
equations for the functions x{z) and y{z): 



dx _ ^l + {d x/dzY 
dz 



= 7T^ ' (41) 



dy 1 d j y'^-idxldz) , , ^. 

y-^'-n. = T ■ — I / „ I , (42) 

Jl + {dx/dz) 



dz a dz 



where z = - 0{t* - t)-^ . 

However, with this choice of parametrization of the curve, the obviously existing sym- 
metry of the system ( |5^ with respect to rotation in the x-z plane is hidden. For taking 
advantage of this symmetry, cylindrical coordinates are more appropriate, with the angle 
coordinate ip serving as the longitudinal parameter: 

{X,Y,Z) = {R{(f, t) cos !f,Y {if, t),-R{(f,t) sin (43) 
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Instead of the equations of motion (p6|-^), we obtain the equivalent system (where a same 
time rescahng as in (|36| - |37D is performed) 



(2 - a)RR 
-{2-a)RY 



1 / Y"R' 



1 RY' 



(44) 
(45) 



Here (..)' = 9. 



This system follows from the Lagrangian written in cylindrical coordinates 



r f Y^ \ 

£^~y (2-«)— F - —^R^^RjAd^. 



The self-similar substitution 

i?(V9, t) = {t* - tfr{^), Y{^, t) = {t* - tfy{v) 



(46) 



(47) 



does not change the meaning of the angle coordinate (p. It leads us to the following pair of 
equations for the functions r{ip) and yi^p): 



,/2 



,1-a 



1 / y°'r' 



ry 



J2 
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(48) 
(49) 



We observe that there is no explicit dependence on in these equations. This property 
helps us to integrate the system. The general solution can be represented in the following 
parametric form (see Appendix B for a detailed derivation): 



(/9(p) = 9?o + arctan(j9) — 



a{\ — a) 



\ (2 -«)(! + «; 



■ arctan p 



a(2 — a) 



y{v) = c- 



r{p) = C2-°= 



Va(2 - a) 



2{2-a) 



a-l 
'2(2-a) 



P 



1 +p2 



(50) 



(51) 



(52) 



a{2 — a) 

where the parameter p runs between the limits — oo < p < +oo, C and (po are arbitrary 
integration constants . The constant C determines the asymptotic slope of the curve at 
large distances from the origin: y k, r/C when r — »• cxo, while the constant (p^ reflects the 
mentioned symmetry of the system with respect to rotations in x-z plane. The condition 
(p5|) for applicability of the local approximation (p^ is satisfied if C^fa ^1. A typical 
self- similar solution x(z) is shown in Fig.^. 
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Figure 3: Self-similar solution x{z) for C = 50, a = 0.1. 



It is interesting to note that the total angle A(f between two asymptotic directions in 
x-z plane does not depend on the parameter C in the long-scale local approximation used 
above: 



A^J = TT 1 



a(l — a) 



\ (2-«)(l + a) 



(53) 



At small values of a, this angle approaches vr. Another remark about A(p is that the above 
expression assumes identical values at a and at 1 — a, so the value (5 = 1/2 results in the 
extremum Acpmin = 27r/3. For this case, the curve lies on the cone y = r/C. 



5 Discussion 

We observed that in the systems ( p4D with < a < 1, finite time singularity formation 
is possible in the self-similar regime. Inasmuch as the condition (^) for applicability of 
the approximate Hamiltonian (153) is satisfied in a range of the parameter C related to the 



self-similar solutions (|50| - |52D , we conclude that in the systems (|^) self-similar collapse of 
two symmetrical singular vortex filaments can also take place. The principal question is 
whether this is also possible for filaments having finite width. If yes, then such solutions are 
analogous to the assumed self-similar solutions of the Euler equation. Though the exponent 
jS ( pO]) differs from 1/2, the difference is small if a is small. However, an important difference 
exists between infinitely thin filaments and filaments with finite width: inside the latter, 
longitudinal fiows take place, caused by a twist of the vortex lines constituting the filament. 
Those fiows keep the width homogeneous along the filament if a local stretching is not 
sufficiently fast. This mechanism acts against singularity formation and, probably, in some 
cases it can prevent a singularity at all. [It is worth mentioning here that for finite width 
vortex structures in the Navier-Stokes equation frame, the usual "outcome" result of the 
Crow instability is vortex line reconnection [^.] Thus, a more-less consistent analysis of 
the general situation should take into account, besides the dynamics of a mean shape of 
the filament, at least the dynamics of the width and the conjugated dynamics of the twist. 
Clearly, we do not need to consider a ^ systems, when we deal with non-singular vortex 
filaments. It should be emphasized that an attempt to take account of finite width of the 
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filament by simply using regularized Green's functions such as G^{r) ~ l/^/rM-~? with 
a constant e, giving correct results for long scale limit of the linearized problem, fails to 
describe the dynamics in the highly nonlinear regime. 

Also, we would like to note that a local approximation analogous to (|3^) is possible for 
arbitrary Green's function GmI'^)- The corresponding long scale Hamiltonian has the form 



where the (positive) function Fm{Y) is related to the function G'Af(r) in the following way: 



The stationary motion with a constant coordinate Yq = b/2 is unstable if the second deriva- 
tive of the function Fm is negative at that value: F^(6/2) < 0. We believe that such 
systems can exhibit locally self-similar collapse, if the asymptotics of the function Fm{Y) is 
power-like at small Y: Fm ~ Y"', with < a < 1. 

The final remark concerns the possibility of including effects caused by the fact that the 
unstable range is finite in the wave-number space into the approximate long scale theory. 
Especially this is important for the case of small values of a, because in that limit the 
condition (^) for applicability of the Hamiltonian ( [3^ ) becomes too restrictive. The idea 
how to improve the approximation is the following. In general, the exact expression for the 
Hamiltonian of a pair of singular filaments, after integration by parts, can be represented 
as the half of the integral over a surface S drawn between the filaments (one half since we 
consider only one from two symmetric strings): 



because the canonical momentum field p created by filaments is determined by a multi- 
valued scalar potential <l>(r): p = V$, which has the additive increment F = /(p • dl) after 
passing around a filament. Also the equality div v = is important for derivation of the 
last expression. In the case of small a, we should just more carefully take account of the 
contribution to the surface integral from the vicinity of filaments. As the result of such 
consideration, we find that for a better approximation it is sufficient to replace in (|3^) the 
projection of the arc-length element by the entire arc-length element and, correspondingly, 
use the Hamiltonian 



We stress once more here that this expression is valid only in the case a ^ 1, F/L ^ 1. 
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Fm{Y) = JJ^ (Gm(0 - GM{^Je + {2Yy)) 





(54) 
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Appendix A 

In order to have some closed expressions for the functions Aa{k) and Ba{k) instead of the 
integral representations (123) and (E^), let us use the following mathematical relations PS[: 



(C + i)^ " 2 r( , 



J(l) 



£"^cosikbC)^ = (6fc)-"cos (^) r(a), (56) 



r+°° cos(qC)dC a/tt /g\^"5 ^ ^ ^ ^ 

where T{x) is the Gamma function, K^{x) is the modified Bessel function of the second kind. 
The integral (|58D results in the equahties 

- f- cos(A:&C)ciC _ A ("^ . ^59^ 

-Jo (e + i)^-r(if^) U; ^''^ 



(C^ + 1)^ r(^)V2 

^cos{kb()d( ^ I kb 



2- 



J") = r""^°^m)dc ^ g (61) 



/o (e + ijT" r 

Thus, we have from (|23) and (ESl): 



5— a 
2 



A„(A;) = 2(1 - a)26"-2j(3) _ 2k%''J^^^ + 2(1 - a)^"-^ (/3_, - , (62) 
B^{k) = 2(1 - a)26-2/(3) _ 2(1 - a)(3 - a)6"-=^ (2(j(5) + 1,.^) - Z^) _ jg,^) . (63) 

Appendix B 



In this Appendix, we explain how the solution (pU|-^) of the system (^^-^) can be obtained. 
Let us introduce the designations 

Q = {dr/dipf, s = r^, (64) 

then consider temporary s as independent variable, and rewrite Eq.(ED[) as follows: 
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or equivalently 



y=-\a—y Js + Q-—' 



ds \^/s + Q 



a \ ds' 

Substituting into this equation the relation 

Q = ,(,^2(1-.) _ ^) 

which follows from the equation (^H]), we have the following equation for y{s): 

a ds ^ 



(66) 



(67) 



(68) 



This first order differential equation is linear for the inverse dependence s{y), and its general 
solution is 



^« , (69) 

a[2 — a) 

where C is an arbitrary constant of integration. Thus, we have the relation between y and 
s = r^. To obtain another relation, between y and if, let us use the equation 



dip 



ds 



(70) 



which gives us the integral 

s'{y)dy 



2(2-a) _(_ (l-2a)(l-a) ^| dy 



a{2-a) J y 



2s{y)^sWi^--)-l J (C^y2i2-.) - ^C-y2i2-.) - ^ 



arctan 



Q2y2{2-a) 



(1 -a2) 



a(l — a) 



a{2-a)J \\{2- a){l + a) 
After introducing the new parameter 



arctan , 



'a{2 — a] 



(J2y2{2^a) _ I 



P 



\ 



C2y2{2-a) _ 11 ^ 

, a(2-a), 



(71) 



(72) 



we arrive at solution of the system (^8|-^) in the form ( ^0| - |5^ . 
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